Non-conformality and non-perfectness of fluid near phase transition 
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We investigate the thermodynamic and transport properties of the real scalar field theory at 
weak as well as strong couplings in the Hartree approximation of Cornwall- Jackiw-Tomboulis (CJT) 
formalism. To our surprise, we find that near phase transition, all the thermodynamic and transport 
properties of the simplest real scalar model at certain strong coupling agree well with the lattice 
results of the complex QCD system. We also demonstrate that the system near phase transition is 
non-conformal, and behaves as a non-perfect fluid. 
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Studying Quantum chroniodynamics (QCD) phase 
transition and properties of hot quark matter at high 
temperature has been the main target of heavy ion col- 
hsion experiments at the Relativistic Heavy Ion col- 
lider (RHIC) and the forthcoming Large Hadron Colhdcr 
(LHC). Before RHIC was turned on, it was expected 
that deconfined quark matter should behave like a gas 
of weakly interacting quark-gluon plasma (wQGP). The 
perturbative QCD calculation gives a large shear vis- 
cosity in the wQGP with 77/s ~ 0.8 for a^ = 0.3 [l|]. 
Surprisingly, in order to fit the elliptic flow at RHIC, the 
hydrodynamic simulation shows that a very small shear 
viscosity is required [4]. Since then it has been believed 
that the system created at RHIC is a strongly coupled 
quark-gluon plasma (sQGP) and behaves like a nearly 
"perfect" fluid [l,i]. 

In fluid dynamics, there is another important trans- 
port coefhcient, the bulk viscosity C, which has often 
been neglected in hydrodynamic simulation of nuclear 
collisions. The zero bulk viscosity is for a conformal 
equation of state and also a reasonable approximation 
for the weakly interacting gas of quarks and gluons. 
For example, the perturbative QCD calculation gives 
C/s = 0.02^2 fo^ 0.06 <as< 0.3 Q. 

Lattice QCD calculation confirmed that 77/s for the 
purely gluonic plasma is rather small and in the range 
of 0.1 — 0.2 Q. However, recent lattice QCD results 
showed that the bulk viscosity over entropy density ratio 
^/s rises dramatically up to the order of 1.0 near the 
critical temperature Tc [7|, |8| . The sharp peak of bulk 
viscosity at Tc has also been observed in the linear sigma 
model Q, and the increasing tendency of C/s below Tc 
has been shown in a massless pion gas [l^l- The large 
bulk viscosity near phase transition is related to the non- 
conformal equation of state [ll|, \1M ■ 

Due to the complexity of QCD in the regime of strong 
coupling, results on hot quark matter from lattice cal- 
culation and hydrodynamic simulation are still lack of 
analytic understanding. In recent years, the anti-dc Sit- 
ter/conformal field theory (AdS/CFT) correspondence 
has generated enormous interest in using thermal A/" = 4 
super- Yang-Mills theory (SYM) to understand sQGP. 
The shear viscosity to entropy density ratio 77/s is as 



small as l/47r in the strongly coupled SYM plasma 
[l3l . [l4| . However, a conspicuous shortcoming of this 
approach is the conformality of SYM: the square of the 
speed of sound c^ always equals to 1/3 and the bulk vis- 
cosity is always zero at all temperatures in this theory. 
Though C/s at Tc is non-zero for a class of black hole 
solutions resembling the e quat ion of state of QCD, the 
magnitude is less than 0.1 [l5[, which is too small com- 
paring with lattice QCD results. 

It has been found in Ref. [1^1 that in the simplest 
real scalar model with Z{2) symmetry breaking in the 
vacuum, 77/s behaves the same way as that in systems of 
water, helium and nitrogen in first-, second-order phase 
transitions and crossover [17| . In this letter, we investi- 
gate the equation of state and bulk viscosity in the real 
scalar model, and compare the result in this simplest rel- 
ativistic system with that of the complex QCD system. 

The Lagrangian of the real scalar field theory has the 
form of 



^ = l^idf.^ 



r^ 



(1) 



with a the mass square term and b the interaction 
strength. This theory is invariant under (p ^ —(p and 
has a Z2 symmetry. In the case of a < and 6 > 0, 
the vacuum at T = breaks the Z2 symmetry sponta- 
neously. The Z{2) symmetry will be restored at finite 
temperature with a second-order phase transition. 

At finite temperature, the naive perturbative expan- 
sion in powers of the coupling constant breaks down. A 
convenient resummation method is provided by the ex- 
tension of Cornwall- Jackiw-Tomboulis (CJT) formalism 
|18| to finite temperature. The CJT formalism is equiva- 
lent to the (f>-functional approach of Luttinger and Ward 
|19| and Baym [20| ■ In our calculation, we only perform 
the Hartree approximation for the effective potential, i.e, 
only resum tadpole diagrams self-consistently and ne- 
glect the exchange diagrams. The effective potential in 
the CJT formalism reads |2li 
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^[<I>.S] = ^J [\nS-HK) + S-\K)SiK)-l] 

+ V2[^,S] + Ui^) , (2) 



where U{4>) — a/2 cj)^ + 6/4 0^ is the tree-level potential, 

and the 2PI potential V2[4>,S] = f 5 (/^ SiK,^))"^ in 
the Hartrce approximation. S{So) is the full(trcc-lcvcl) 
propagator and takes the form of S~^{K,(f>) = —K^ + 



m2(0), So'{K, 



-K + rnQ{4>) with the tree- level 



mass TOq ~ a + 3b (j) . 

The gap equations for the condensation (Jjq and scalar 
mass m are determined by the self-consistent one- and 
two-point Green's functions 
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The entropy density is determined by taking the 
derivative of effective potential with respect to temper- 
ature, i.e, s = —d^{(j)t)) / dT . In the symmetry break- 
ing case, the vacuum effective potential or the vacuum 
energy density is negative, i.e, Q.^ = Q,((f)Q)\T=a < 0. 
As the standard treatment in lattice calculation, we in- 
troduce the normalized pressure density pr and energy 
density er as pr = —ftr and ex = —pr + Ts, with 
VIt = f^((/'o) — r^ti.The equation of state pt{^t) is an 
important input into hydrodynamics. The square of the 
speed of sound c^ is related to pr /^t and has the form 
of 



dp s s 

d^ ~ Tds/dT ~ a ' 



(4) 



where Cy = de/dT is the specific heat. At the criti- 
cal temperature, the entropy density as well as energy 
density change most fastly with temperature, thus one 
expect that c^ should have a minimum at Tc. The trace 
anomaly of the energy-momentum tensor T^'^ 



r^'^' _ er- ipr 
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(5) 



is a dimensionless quantity. We define A/d as the "in- 
teraction measure" , with d the degeneracy factor. 

The bulk viscosity is related to the correlation function 
of the trace of the energy-momentum tensor 6*1^: 



1 1 r°° 

C = - lim - / dt 
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According to the result derived from low energy theorem, 
in the low frequency region, the bulk viscosity takes the 
form of 
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{-16eT + 9T5 + Ta + 16|e„|} . (7) 



with the negative vacuum energy density e-u = fly — 
ri((/)o)|T=0; and the parameter ujq = uJo{T) is a scale at 
which the perturbation theory becomes valid. 

The conformal limit has attracted much attention 
in recent years, since people are trying to under- 
stand strongly interacting quark-gluon plasma by using 




FIG. 1: The interaction measure {er — Spt)/T as a function 
of temperature T/Tc for different coupling strength b. 



AdS/CFT techniques. In conformal field theories includ- 
ing free field theory, pt/^t — c^ = 1/3, A = 0, and the 
bulk viscosity is always zero. Lattice results show that 
at asymptotically high temperature, the hot quark-gluon 
system is close to a conformal and free ideal gas. 

However, lattice results show that near deconfine- 
ment phase transition, the hot quark-gluon system de- 
viates far away from conformality. Both pt/^t and c^ 
show a minimum around 0.07, which is much smaller 
than 1/3. For the SU{S) pure gluon system[ll|, the 
peak value of the trace anomaly A^^^. reads 3 ~ 4 
at Tfnax and the corresponding "interaction measure" 
is A^^y/dc = 0.2 ^ 0.25, with the gluon degeneracy 
factor do = 16. (Note that here Tmax — 1-lTc is the 
temperature corresponding to the sharp peak of A.) For 
the two-fiavor case [12], the lattice result of the peak 



value of the trace anomaly A 



Nf=2 
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reads 8 '^ 11, the 



corresponding interaction measure at Tmax is given as 
^LAT /(^G + dq) = 0.28 ~ 0.4, with quark degeneracy 
factor dq = 12. There have been some efforts trying 
to understand trace anomaly in gluodynamics near and 
above Tc in terms of dimension t wo g luon condensate 
and an effective "fuzzy" bag model [23 |. 

As a reference for complex QCD system, we investi- 
gate the equation of state and transport properties for 
the real scalar field theory with Z(2) symmetry breaking 
in the vacuum and 2nd order phase transition at finite 
temperature. The trace anomaly A the specific heat Cy 
as well as bulk viscosity to entropy density ratio C/s show 
upward cusp at T^ and their peak values increase with 
the increase of coupling strength. The ratio of pressure 
density over energy density pt/^t and the square of the 
sound velocity c^ show downward cusp at Tc, which is 
similar to the behavior of r)/s found in Ref. |16| . and 
the cusp values decrease with the increase of coupling 
strength. These cusp behaviors at phase transition re- 
semble lattice QCD results. In Figs. [Tl[21 and [3] we only 
show the trace anomaly A and the specific heat Cy and 
the ratio of bulk viscosity to entropy density ratio (^/s 
as functions of T/Tc for different coupling strength b. 

In the weak coupling case when b ~ 0.3, the cusp val- 
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llll |12| | , the Polyakov-loop Nambu-Jona-Lasinio (PNJL) 
model HI,!!!, and black hole duals EM- 
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FIG. 2: The specific heat Cv as a function of temperature 
T/Tc for different coupling strength b. 
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FIG. 3: The bulk viscosity over entropy density ratio (^/s as a 
function of temperature T/Tc for different coupling strength 
b. We have used uja = lOT. 



ues of pt/et and c^ at Tc are close to the conformal value 
1/3, both the trace anomaly A and the bulk viscosity 
to entropy density ratio C/s a-t Tc arc close to confor- 
mal value 0. However, the shear viscosity over entropy 
density ratio 77/s is around 2000, which is huge compar- 
ing with the AdS/CFT limit I/Att. Here we have used 
the method in Ref.[16| to derive rj/s. To our surprise, 
we find that when 6 = 30, the strongly coupled scalar 
system can reproduce all thermodynamic and transport 
properties of hot quark-gluon system near Tc- pt/^t at 
Tc is close to the lattice QCD result 0.07, A/d = 0.48 
(d = 1 for scalar system) at Tc is close to the lattice re- 
sult of the peak value A^^^ /{da + dg) ~ 0.4 at T^ax- 
The bulk viscosity to entropy density ratio C/s at Tc is 
around 0.5 ^ 2.0, which agrees well with the lattice re- 
sult in Ref. @]. (Note, here ^/s = 0.5,2 correspond 
to ujo = lOT, 2.5T, respectively.) More surprisingly, the 
shear viscosity over entropy density ratio rj/s at Tc is 
0.146, which also beautifully agrees with lattice result 
0.1 - 0.2 in Ref. @. 

In Table m we compare our results of equation of state 
and transport properties in scalar field theory at Tcjand 
corresponding results in lattice QCD calculations [3, [^, 
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TABLE I: Thermodynamic and transport properties at 
T/Tc = 1 in scalar theory at weak couplin g b — 0.3 and 
strong coupling b = 30, in lattice QCDffl, @,[l3, [12, PNJL 
model [23 . [24J . and black hole dules [15| . The degeneracy 
factor d= \ for real scalar model. 

In summary, in the Hartree approximation of C JT for- 
malism, we have investigated the equation of state and 
transport properties of the real scalar field model with 
Z(2) symmetry breaking in the vacuum and 2nd-order 
phase transition at finite temperature. 

We have seen that at phase transition, the system ei- 
ther in weak coupling or strong coupling shows some 
common properties: 1) pT/tT, the square of the speed 
of sound cf as well as rj/s exhibit downward cusp behav- 
ior at Tc- 2) The trace anomaly A, the specific heat Cy 
as well as C,/ s show upward cusp behavior at Tc- The 
cusp behavior is related to the biggest change rate of 
entropy density at Tc- 

At weak coupling, the scalar system near phase tran- 
sition is asymptotically conformal. However, the shear 
viscosity is huge. At strong coupling, the scalar system 
near phase transition is highly non-conformal, the shear 
viscosity is small, but the bulk viscosity is large. Be- 
cause of the high non-conformality near phase transition, 
AdS/CFT method maybe cannot help us understand the 
strongly interacting quark gluon plasma. Unexpectedly, 
the simplest scalar field model can do the job. We have 
found that lattice QCD results on the equation of state 
and transport properties near phase transition can be 
amazingly very well described by the simplest real scalar 
model at strong coupling when b = 30. 

Small shear viscosity means that the system is strongly 
coupled. At the same time, strongly interacting system 
means large non-conformality and large bulk viscosity. 
Therefore, the strongly interacting system cannot be per- 
fect fluid, especially near phase transition. It is urgent 
to include the bulk viscosity correction and the noncon- 
formal equation of state in hydrodynamics to investigate 
hadronization and freeze-out processes of QGP created 
at heavy ion collisions [2^, [20] , it would be also interest- 
ing to investigate how bulk viscosity affects charm radial 
flow [23] at RHIC. The universal property of large bulk 
viscosity near phase transition r nay play important role 
in the evolution of early universe [2g| , and also may affect 



the cooling of neutron stars [29|. 

At the end, we hope to point out some hniitations of 
our results: Firstly, the results of thermodynamic prop- 
erties in this paper are based on Hartree approxima- 
tion in the CJT formalism. As we know that mean- 
field approximation cannot describe critical phenom- 
ena very well. For 2nd-order phase transition in Z(2) 
model, the specific heat C„ should diverge at the critical 
point, and behave as t~" near the critical point, with 
t = {T — Tc)/T(. and a = 0.11. However, in Hartree ap- 
proximation of CJT formalism, though we observe the 
weak divergence of C^ at T^ in the case of strong cou- 
pling, we can only see a weak upward cusp of C„ at Tc 
in the case of weak coupling. Secondly, the results of 
bulk viscosity in this paper are based on Eq. ^. The 
limitation of Eq. ([7]) has been analyzed in Rcfs. [3C| 
and [3l|. From Eq. ([7]), we see that the bulk viscosity 
is dominated by C^ at Tc- If C^ diverges at Tc, the bulk 
viscosity should also be divergent at the critical point 
and behave as <^". However, the detailed analysis in the 



Ising model in Ref. [32| shows a very different divergent 
behavior C, ~ i^^'^"'"", with z ~ 3 the dynamic critical ex- 
ponent and v ~ 0.630 the critical exponent in the Ising 
system. Thirdly, we should keep in mind that our results 
at strong coupling are from effective theory. From renor- 
malization analysis, the scalar theory will hit a Landau 
pole when jSb > b with the (3 function (3h — 96^/167r^. 
The results for large h in the scalar theory are not guar- 
anteed to be valid in the CJT formalism. 
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